February 5, 2008 3:48 Proceedings Trim Size: 9.75in x 6.5in 



procfachruddin 



THREE-NUCLEON SCATTERING AT INTERMEDIATE ENERGIES 



I. FACHRUDDIN 

Departemen Fisika, Universitas Indonesia, Depok 16424, Indonesia 

CH. ELSTER 

Department of Physics and Astronomy, Ohio University, Athens, OH 34701, USA 

W. GLOCKLE 

Institut fur Theoretische Physik II, Ruhr-Universitdt Bochum, D-44780 Bochum, 

Germany 



By means of a technique, which does not employ partial wave (PW) decompositions, 
the nucleon-deuteron break-up process is evaluated in the Faddeev scheme, where only 
the leading order term of the amplitude is considered. This technique is then applied 
to calculate the semi-exclusive proton-deuteron break-up reaction d(p, n)pp for proton 
laboratory energies Ei a ^ of a few hundred MeV. A comparison with PW calculations 
is performed at 197 MeV projectile energy. At the same energy rescattering processes, 
which are not included in the 3D calculations yet, are shown to be still important in 
the full Faddeev PW calculations, especially for the cross section and the analyzing 
power A y . Next kinematical relativistic effects are investigated for projectile energies 
up to about 500 MeV. At the higher energies, those relativistic effects start not to be 
negligible, especially in the peak of the cross section. 



1. Introduction 

In order to investigate the short distance behavior of three nucleon (3N) as well 
as nucleon-nucleon (NN) forces it is necessary to consider 3N scattering at higher 
energies of a few hundred MeV. The partial wave (PW) technique, which takes 
as basis states a specific number of angular momentum eigenstates, has a long 
history of being employed to solve Faddeev equations for 3N scattering [T] . However 
as the energy increases the number of contributing angular momenta proliferates, 
leading to increased algorithmic difficulties as well as more tedious algebraic work. 
Therefore, a new and alternative approach is needed, which is not based on PW 
basis. 

In momentum space the choice is to work directly with momentum vector states. 
For three boson scattering this approach has been pioneered successfully carried out 
in Refs. [2]. We develop our approach [3], referred here as the 3D technique, for 
three nucleon scattering, allowing the use of realistic nuclear forces, since spin and 
isospin are also taken into account. Since the NN system serves as input for our 
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3N calculations, we successfully applied the 3D technique to both NN scattering 4 
and the deuteron [5] . 

We evaluate the nucleon-deuteron (Nd) break-up process using the Faddeev 
scheme. However, we do not solve the full Faddeev equations but rather concentrate 
on the first order term in the multiple scattering series given by them. Thus, we 
assume that at intermediate energies of a few hundred MeV the leading term may 
sufficiently describe the scattering process. We include relativistic effects in the 
kinematics to investigate their size as function of energy. For our application we 
concentrate on the semi-exclusive d(p, n)pp reaction and calculate the spin averaged 
differential cross section, the neutron polarization, the proton analyzing power and 
the polarization transfer coefficients. There are in fact experimental data in the 
energy range up to about 500 MeV [5]-[S]. The 3D technique works with any 
potentials given in operator form, not the PW projected ones. We employ the NN 
potentials AV18 [9J and Bonn-B [TU]. 

2. The Nd Break-Up Amplitude 

The break-up process is given by the Nd break-up operator Uq given by 

U = (1 + P)TP, (1) 

where P = P12P23 + P13P23 is a permutation operator, T stands for NN t-matrix. 
The term TP is the leading term of the 3N transition operator Tp obeying Faddeev 
equation Tp = TP + TGoPTp, with Go being the free three-nucleon propagator. 
The operator Uq is fully anti-symmetrized and can be written as a sum of three 
terms in TP, the first of which is £/ (1) = TP. The other terms f/ (2) = P12P23TP 
and Uq 3 " 1 — P12P23TP are related to Uq^ by means of permutations. 

Without employing PW decompositions, except for the deuteron, the Nd break- 
up amplitude Uo(p, q) is given in Eq. (J2j) (see Ref. [3] for derivation), where m s i, Ti 
(z = 1,2,3) are the spins and isospins of the three nucleons in final state, mP sl , t® the 
spin and isospin of nucleon 1 acting as the projectile, ^ d d the deuteron state, M<j 
the projection of total angular momentum of the deuteron along an arbitrary z-axis, 
tpi(7T r ) the deuteron wave function components (s and d waves), T^fJ{p, ir, cos 9'; E p ) 
the NN t-matrix elements in a momentum-helicity basis [4] for given parity 7r, total 
spin S and isospin t, final and initial helicities A and A', and a center of mass 
energy E p of the 23-subsystem, d A , A (9)'s being the rotation matrices [11] . The 
Jacobi momenta p and q describe the 3N kinematics in the final state, where p 
is the relative momentum between nucleon 2 and 3, q the relative momentum of 
nucleon f to the 23-subsystem, and qo the relative momentum of the projectile to 
the deuteron. 

The operator C/o(p, q) given in Eq. ([2]) is derived within the framework of the 
nonrelativistic Faddeev scheme. By adopting the formulation given in Ref. [12] . 
relativistic kinematics is introduced. Thus, we re-evaluate the Jacobi momenta, 
carry out corresponding Lorentz transformations to the two- and three-particle cm. 
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subsystems, and employ relativistic energy-momentum relations. To calculate the 
observables, a relativistic description of the cross section is employed. 



^o(p,q) = ( pqm s im S 2m s3 TxT2T3 



( ) 2+Tl j _ „-irA <6„-AU,")^t f l 1 
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+ ^ 2) (p,q) + C/f (p,q), (2) 



with 



1 1 
7r=2 q + q ° ^r^-q-^^o (3) 

sin 6 p sin ft^ cas(<f> p — (j)^) (4) 

(5) 



EN=- S e iN ^-^d s NA (e p )d s NA ,(e^ 



„j(AV'-An) _ Z^Af=-S e "WAV^P 



3. The Semi-exclusive Proton-Deuteron Break-Up Reaction 

Our formulation is applied to the semi-exclusive proton-deuteron (pd) break-up 
reaction d(p, n)pp. We calculate the spin averaged differential cross section, the 
neutron polarization, the proton analyzing power and the polarization transfer co- 
efficients for proton laboratory energies Ei ab up to about 500 MeV. 

First, we perform comparisons with calculations based on the well established 
PW technique, which also include only the leading term of the Faddeev amplitude. 
Both schemes agree for projectile energies Ei^ < 200 MeV. However, at about 200 
MeV deviations occur in the cross section peak, as shown in Fig. [T] for the cross 
section at Ei a b = 197 MeV, where the PW calculation has not sufficiently converged 
to the 3D calculation. The PW calculation shown in Fig. [1] takes into account NN 
angular momenta up to j = 5 and 7, and 3N states for total 3N angular momenta 
up to J — 31/2, which is a technical maximum at present. 

Since we do not take rescattering terms into account, we nevertheless have to 
estimate their effects. At E^ = 197 MeV we do this by comparing the 3D calcula- 
tion with a full Faddeev PW calculation, and display the results in Fig. O We see 
that the higher order multiple scattering contributions lower the cross section and 
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Figure 1 . The spin averaged differential cross 
section at E lah = 197 MeV, 9 = 13°, based on 
the Bonn-B potential. 
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are essential for the description of the 
analyzing power especially for the small 
energies of the outgoing neutron. 

Next, we investigate effects of rela- 
tivistic kinematics on semi-inclusive ob- 
servables as function of the projectile 
energy. In Fig. [3] we compare our 
nonrelativistic 3D calculations with the 
corresponding relativistic ones, and see 
that the position as well as the height 
of the quasi-free peak are influenced. 
Since the position of the quasi-free peak 
is solely determined by kinematics, it is 
satisfying to see that the use of rela- 
tivistic kinematics put the calculated 
peak at the right position with respect 
to the data. As expected the effect 
of relativistic kinematics increases with 
increasing projectile energy. 



0.16 



«qf -0.07 - 




-0.30 1 1 1 

70 110 150 

E n [MeV] 



Figure 2. The spin averaged differential cross section (a) and the analyzing power A y (b) at 
Elab = 197 MeV, 8 = 37°. The calculations are based on the AV18 potential, the data from 
Ref. [8]. 
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4. Summary and Conclusions 

We calculated the semi-exclusive pd break-up process within the Faddeev scheme 
in the leading order in a multiple scattering expansion in a 3D formulation based 
directly on momentum vectors. This technique has proved to be a viable alternative 
to traditional PW decompositions. For energies smaller than 200 MeV our 3D 
calculations show perfect agreement with PW calculations. At energies larger than 
200 MeV the inclusion of relativistic kinematics proves essential to obtain the correct 
position of the quasi-free peak in the (p,n) charge exchange reaction. 
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